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ABSTRACT 

We have observed a sample of 35 long-period variables and four Cepheid variables in the vicinity of 
23 Galactic globular clusters using the Infrared Spectrograph on the Spitzer Space Telescope. The long- 
period variables in the sample cover a range of metallicities from near solar to about l/40th solar. 
The dust mass-loss rate from the stars increases with pulsation period and bolometric luminosity. 
Higher mass-loss rates are associated with greater contributions from silicate grains. The dust mass- 
loss rate also depends on metallicity. The dependence is most clear when segregating the sample by 
dust composition, less clear when segregating by bolometric magnitude, and absent when segregating 
by period. The spectra are rich in solid-state and molecular features. Emission from alumina dust 
is apparent across the range of metallicities. Spectra with a 13-/^m dust emission feature, as well as 
an associated feature at 20 /iin, also appear at most metallicities. Molecular features in the spectra 
include H2O bands at 6.4-6.8 /mi, seen in both emission and absorption, SO2 absorption at 7.3-7.5 fim, 
and narrow emission bands from CO2 from 13.5 to 16.8 /jja. The star Lynga 7 VI has an infrared 
spectrum revealing it to be a carbon star, adding to the small number of carbon stars associated with 
Galactic globular clusters. 

Subject headings: globular clusters: general — stars: AGB and post-AGB — infrared: stars — cir- 
cumstellar matter 



1. INTRODUCTION 

Stars ascending the asymptotic giant branch (AGB) 
burn hydr ogen and helium in sh ells around an inert C-0 
core (e.g., Ilben fc Renz ini 1983). The fusion of helium 
to ca rbon proceeds by the triple-o sequence (|Salpeterl 
1952) in thermal pulses which lead to the dredge-up of 
freshly produced carbon to the surface of the star. If the 
envelope of the star is not too massive and the dredge- 
ups are sufficiently strong, then enough carbon reaches 
the surface to drive the photospheric C/O ratio over 
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unity. Generally, the formation of CO molecules will 
exhaust whichever of the carbon or oxygen is less abun- 
dant, leading to a chemical dichotomy. Carbon stars 
produce carbon-rich dust, and oxygen-rich stars produce 
oxygen-rich dust. Stars in the range from ^2-5 Mq be- 
come carbon stars, although both the lower and upper 
mass limits decrease in more m etal-poor environments 
(e.g. iKarakas fc Lattanzioll2007fh 

Several su r veys w ith the Infrared Spectrograph (IRS; 
iHouck et alj 2004) on the Spitzer Space Telescope 
(jWerner et al.ll2004l ) have probed how the production of 
dust by evolved stars depends on metallicity by observ- 
ing AGB stars and supergiants in the Magallanic Clouds 
and other nearby Local Group galaxies. Carbon stars 
dominate these samples. The rate at which they pro- 
duce dust does not vary significantly with metallicity 
(jGroenewegen et alj I2007t ISloan et ail 120081 and refer- 
ences therein). It would appear that they produce and 
dredge up all of the carbon they ne ed to form dust, re- 
gardless of their initial abundances ()Sloan et alj |2009). 

The published Local Group samples contain fewer 
oxygen-rich evolved stars, maki ng any conclusions about 
this population more tentative. ISloan et al.l (|2008h com- 
pared oxygen-rich sources in the Galaxy, Large Magel- 
lanic Cloud (LMC), and Small Magellanic Cloud (SMC), 
and they found that as the metallicity of the sample 
decreased, the fraction of stars with a dust excess de- 
creased. Infrared photometric surveys of several dwarf 
irregular galaxies in the Local Gr oup also show a tr end 
of less dust at lower metallicities (jBover et a l. 2009a) . 

The spectroscopic evidence for a dependency of dust 
production on metallicity is not strong. Stars in the 
Magellanic Clouds can pulsate with periods of ^700 days 
or longer, and these longer-period variables are usually 
embedded in significant amounts of dust, with no ap- 
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parent dependence on metallicity. Additional oxygen- 
rich sources from the SMC fu rther blur the trends no- 
ticed before (jSloan eFaLl[2010). One major problem with 
these samples is that oxygen-rich evolved stars can come 
from three distinct populations. The red supergiants and 
intermediate-mass AGB stars (or super-AGB stars) are 
too massive to become carbon stars, while the low-mass 
AGB stars have too little mass. This mixing of pop- 
ulations confuses the observed samples of oxygen-rich 
evolved stars in the Galaxy, the Magellanic Clouds, and 
more distant irregulars in the Local Group. 

Globular clusters provide another means of investi- 
gating dependencies of dust production on metallicity. 
Like the Magellanic Clouds, most globular clusters have 
known metallicities, and they are at known distances, al- 
lowing us to directly determine their luminosities. Unlike 
the Magellanic Clouds, globular clusters have old popu- 
lations of stars, with few significantly younger than 10 
billion years. This limits the sample to masses of ~1M 
o r less, below the lower mass limit for carbon stars. 

iLebzelter et al.l ((2006) identified and observed 11 long- 
period variables (LPVs) in the globular cluster 47 Tuc 
with the IRS on Spitzer. Their data were consistent 
with a shift from relatively dust-free AGB stars to more 
deeply embedded sources with higher mass-loss rates at 
a luminosity of ~2000 L . But with a sample of only 
one globular cluster, they were unable to address the im- 
portant question of how the dust properties depend on 
metallicity. 

We have used the IRS on Spitzer to observe a sample of 
39 variable stars in 23 globular clusters spanning a range 
of metallicity from nearly solar to only a few percent 
of solar. This paper presents an overview of the pro- 
gram. §[2] describes the sample of clusters and individual 
stars and the observations, including near-infrared pho- 
tometric monitoring prior to the Spitzer observations, the 
IRS spectroscopy, and nearly simultaneous near-infrared 
photometry. § [3] assesses the membership of our targets 
within the clusters, while § [4] describes the analysis and 
the results. §[5] treats some unusual objects individually, 
and § |6] examines the larger picture of evolution on the 
AGB. 

2. OBSERVATIONS 

2.1. The Clusters 

Table [TJ presents the metallicity, reddening, and dis- 
tance modulus to the clusters in our sample. The values 
for each are based on a review of the literature, start- 
ing with the catalog of lHarrisI (|2003l ). which was pub- 
lished in 1996 and updated on the web in 2003. Most 
of the estimates for metallicity and reddening are un- 
changed between the two editions for our target clusters, 
but most of the distances have changed, usually by small 
amounts. We have also included measurements of the 
metallicity and distance modulus published since 1996, 
taking care not to overweight references included in the 
2003 updates to the bibliography of the catalog. We 
have averaged the measurements for metallicity, redden- 
ing, and distance modulus. For the latter two, we also 
quote the standard deviation as the uncertainty, which 
propagates through the quantities we derive. We take 
the minimum uncertainty for E(B — V) to be 10% of the 
mean. In some cases, we have quoted the uncertainties 



given by the referenced papers. Some data which vary 
substantially from the others have been excluded (and 
are not included in the references in Table [TJ). 

For IC 1276 and Terzan 12, we excluded distance mod- 
ulus measurements inconsistent with our estimates based 
on the period-K relation (§ I3.1|) . For Terzan 5, we report 
the distance modulus based on the analysis of the four 
relatively unreddened variables in our sample (§ 13. ip . 

The metallicity of Palomar 6 is problematic, with mea- 
surements in the refe reed liter ature ranging from —1.08 
(|Lee fc Carnevil2002D to +0.2 (IMinnitil 119951) . Interme- 
diate measurements of —0.74 (jZinnl Il985f ) and —0.52 
( Step hens fc Fr ogcl 2004]) have also been reported. Aver- 
aging these metallicities gives <[Fe/H]> = — 0.54±0.54. 
With a standard deviation as large as the mean, we are 
unable to determine a metallicity for this cluster with 
any confidence. 

The metallicity of the cluster is important, because we 
assign individual stars to bins based on their metallicity. 
Each bin will contain between five and eight targets, all 
with a similar metallicity. By comparing how the prop- 
erties of the stars vary from one bin to the next, we can 
assess how the mass-loss and dust production depend on 
metallicity. 

2.2. Photometric Monitoring and Sample Selection 

We selected the variable stars in the sample based 
on JHK photometry from the South African Astronomi- 
cal Observatory (SAAO). Observations were made using 
the SIRIUS near-infrared camera (Si multaneous-Color 
Infra red Imager for Unbiased Surveys; iNagavama et al.1 
120031 ) at the 1.4-m Infrared Survey Facility (IRSF) 
telescope over the period from 2002 to 2005. The 
IRSF/SIRIUS survey imaged fields in the vicinity of glob- 
ular clusters at J, H, and K s (central wavelengths of 1.25, 
1.63, and 2. 14 /mi) with a typ ical sampling interval of 
40-60 days (|Matsunagall2007al) . 

The field of view of the survey was 7.7 arcmin 2 , and 
most clusters were imaged out to the half-mass radius. 
Many new variables were discovered, including eight 
new long-period variables: Lynga 7 VI, Palomar 6 VI, 
Terzan 5 V5-V9, and Terzan 12 VI. All eight of these 
stars are Miras, and their periods range from 260 to 570 
days. The remaining 31 targets in our sampl e were pre- 
viously known variables (|Clement et al.ll2001l and refer- 
ences therein). 

We determined periods, mean magnitudes, amplitudes 
and phases for each variable by performing Fourier anal- 
ysis and minimum-^ 2 fits to the individual photometric 
observations in each SIRIUS filter for each target. The 
results appear in Tableland Figure [TJ The quoted pe- 
riod and phase are the average from the three filters, 
except in cases where one filter disagreed substantially 
from the other two, or for Palomar 6 VI, which was too 
reddened for measurements at J. Generally, the uncer- 
tainty in the period was less than one davFl We lack 
monitoring observations of NGC 6171 VI, NGC 6356 V5, 
and Palomar 10 V2. For thes e sources, Tabl e JH co ntains 
variability information fr om lOement et al.l (|2001[ ) and 
photometry from 2MASS ([Skrutskie et al.ll2006D . 

The four "CW" variables are Pop. II Cepheids (of the 

13 In three cases, Lynga 7 VI, Terzan 5 V9, and NGC 6838 VI, 
it was larger, 5-6 days. 
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TABLE 1 

Globular clusters in the study 



Globular 
Cluster 


Alternative 
Name 


[Fe/H] 


Ref. a 




E(B- 


-v) 


Ref. a 




Distance 
Modulus 


Ref. a 


NGC 362 






— 1. 


.20 


1,2,3,4,5,6 


0. 


.04 


± 





.02 


1,3,4,7 


14 


.83 


± 


0. 


15 


1,3,4,6,7,8 


NGC 5139 




Cen 


-1. 


.63 


1,5 


0. 


12 


± 





,01 b 


1 


13 


.70 


± 


0. 


17 


1,9,10,11,12,13,14,15,16 


INGC 5904 


M 


5 


-1. 


20 


2,3,4,5,6 


0. 


.03 


± 





.01 


1,3,4,7,17 


14. 


.46 


± 


0. 


10 


1,3,4,6,7,8,17 


NGC 5927 






— 0. 


.35 


1,3,4,5,6 


0. 


.46 


± 


o 


05 b 


1,3 


14. 


.45 


± 


0. 


07 


1,3,4,6 


Lynga 7 






-0. 


66 


1,6,18 


0. 


.73 


± 





.12 


1,18 


14. 


.31 


± 


0. 


10 


1,18 


NPP R1 71 


M 


107 


-0. 


.98 


1 Q A K p. 


n 

u. 


.38 


± 





.08 


1 Q A 1 7 


10 


.89 


± 


0. 


18 


J. ,0,1, 1 / 


"WPP 


M 


10 


— 1. 


51 


1 Q K ft 
1,0,0,0 


u. 


.28 


± 





,03 b 


1.0 


1 Q 
1.1. 


.30 


± 


0. 


12 


1,0 


NGC 6352 






-0. 


.69 


1,3,5,6 


0. 


21 


± 





02 b 


1.3 


13 


.85 


± 


0. 


06 


1,3,17 


NGC 6356 






-0. 


.50 


1,4,5 


0. 


.28 


± 





03 b 


1 


15. 


.93 


± 


0. 


09 


1,4 


NGC 6388 






-0. 


57 


1,4,5,19 


0. 


35 


± 





.04 


1,20 


15 


.47 


± 


0. 


12 


4,20,21 


Palomar 6 










22 


1. 


.36 


± 





,14 b 


1,23,24,25 


14. 


.23 


± 


0. 


10 


1,23,24,25 


Terzan 5 






-0. 


08 


1,5,23,26,27 


2. 


.31 


± 





,23 b 


1,23,27 


14 


.11 


± 


0. 


13 


28 


NGC 6441 






-0. 


.56 


1,4,5,29 


0. 


.50 


± 





,05 b 


1,29 


15 


.57 


± 


0. 


16 


"1 A on fin nl 

1,4,29,30,31 


NGC 6553 






-0. 


.28 


1,3,5 


0. 


.72 


± 





.11 


1,3,23 


13 


.63 


± 


0. 


22 


1,3,23,30 


IC 1276 


Palomar 7 


-0. 


69 


1,5 


1. 


12 


± 





,ll b 


1,32 


13 


.66 


± 


0. 


15 b 


1 


Terzan 12 






-0. 


.50 


23 


2. 


.06 


± 





,21 b 


23 


13. 


.38 


± 


0. 


15 b 


1 


NGC 6626 


M 


28 


-1. 


.46 


1,5 


0. 


.40 


± 





,04 b 


1 


13. 


.73 


± 


0. 


15 b 


1 


NGC 6637 


M 


69 


-0. 


.66 


1,3,4,5,6,33 


0. 


.16 


± 





.02 


1,3,33 


14 


.76 


± 


0. 


09 


1,3,4,21,33 


NGC 6712 






-0. 


.97 


1,3,5 


0. 


.43 


± 





.04 


1,3,4 


14 


.24 


± 


0. 


09 


1,3,4 


NGC 6760 






-0. 


.48 


1,4,5 


0. 


.77 


± 





08 b 


1 


14 


.59 


± 


0. 


29 


1,4 


NGC 6779 


M 


56 


-2. 


05 


1,5,6 


0. 


.19 


± 





.08 


1,34 


15 


.23 


± 


0. 


30 


1,34 


Palomar 10 






-0. 


10 


1 


1. 


.66 


± 





,17 b 


1,31 


13 


.86 


± 


0. 


37 


1,32 


NGC 6838 


M 


71 


-0. 


73 


1,2,3,4,5,6 


0. 


26 


± 





,03 b 


1,2,3 


12. 


.99 


± 


0. 


10 


1,3,4,7 



a Ref erences : (11 iHajrlsj lISOOoT); (2) ICarretta fc Gratton] (Ii"997l 1; (31 IFerraro et al.l (IT999T ); ( 41 IRecio-Blanco et al~1 
(2005); (5) Carrctta et al. (2009); (6) Dottcr et al. (2010); (7) Carretta et al. (2000); (8) Gratton et al. (1997); 
(9) Morgan & Dickcrson (2000) ; (10) McNamara (2000); (11) Thompson et al. (2001); (12) Kaluzny et al. (2002); 
(131 ICaputo et al.l (120*0211 ; (141 IDel Principe et ail 120(5(3') ; (15) ICatelanl (120061); (16) Ivan de Ven et al.l ( 1200611; (17) 
Salaris & Weiss (1998); (18) Saraiedini (2004); (19) Carrctta ct al. (2007); (20) Dalcssandro ct al. (2008); (21) Matsunaga 
(2007 13); (22) See S [2Tl (23) Barbuv ct al. (19951); (24) Ortolani ct al. (199 8]); (25) [Te c fc Carney (2002 ); (26) 
IQriglia & Rich (2004); (27) Val cnti et al.l ( 1200711; (28) This work (5|3~H; (29) Valcnti ct al. (2004); (30) Hcitsch lTRichtlerl 
(119991); (31 ) Matsuna ga et al.l l !2009in ^1 lF3arbuv et al.l (I1998bl) ; (33) [Valenti et al.l ( 1200511 ; (34l lIvanov et al.l (120001 ); (35) 
IKaisler et al.l (|1997l ) 



Uncertainty assumed. 

W Vir type). They all have periods which are very short 
compared to the monitoring period, so any irregularities 
in the pulsations can accumulate to smear out the appar- 
ent periodicity in Figure [T] Furthermore, we are unable 
to reproduce the published periods for NGC 6626 V17 
and NGC 6779 V6, which are 92.1 and 90 days, re- 
spectively (jClement et al.ll2001f ). Instead, we find a pe- 
riod of 48.6 days for NGC 6626 V17 and 44.9 days for 
NGC 6779 V6. These appear to be overtones of the pre- 
viously published periods, but we cannot say whether 
that is the result of an actual mode-shift in pulsation or 
sampling resolution of the older data. 



2.3. IRS Spectroscopy 

We observed the sample of variables in Table [5] with 
the IRS on Spitzer, using the Short-Low (SL) module on 
all targets and the Long-Low (LL) module on all but the 
three faintest. The data reduction began with the flat- 
fielded images produced by the S18 version of the data 
pipeline from the Spitzer Science Center (SSC). To sub- 
tract the background from these images, we generally 
used images with the source in the other aperture in SL 
(aperture differences) and images with the source in the 
other nod position in the same aperture in LL (nod dif- 
ferences). In crowded fields, we occasionally had to re- 
sort to whichever background image resulted in the least 
confusion. We then corrected the bad pixels in the dif- 



ferenced images using the imclean IDL packag43- Bad 
pixels included those flagged in the bit-mask images ac- 
companying the data images and those flagged as rogues 
in the campaign rogue masks provided by the SSC. We 
stacked the campaign rogue masks to produce super- 
rogue masks, considering a pixel to be bad if it had been 
identified as such in any two previous campaign rogue 
masks. 

We used two methods to extract spectra from the cor- 
rected images. The older method uses the profile, ridge, 
and extract routines available in the SSC's Spitzer IRS 
Custom Extractor (SPICE). This method produces the 
equivalent of a tap ered-column extraction in SMART 
(Hig don et al. 20040, 

summing the flux within an aper- 
ture which increases in width proportionally with wave- 
length. It is the basis for the spectra in most of the 
IRS papers referenced in the introduction. This globular 
cluster data-set is one of the first to make use of the op- 
timal ext raction algorithm availab le in the new release of 
SMART (jLebouteiller et al.|[2010t) . This algorithm fits a 
super-sampled point-spread function (PSF) to the data, 
improving the signal/noise (S/N) ratio of the data by a 
factor of ~1.8, typically. Optimal extraction also allows 
the extraction of sources with overlapping PSFs, which 

14 Available from the SSC as irsclean. 

15 SMART is the Spectroscopic Modeling, Analysis, and Reduc- 
tion Tool. 
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TABLE 2 

Observed variables in globular clusters 





RA Dec. 


IRS observations 


Var. 


Period 


Mean magnitudes 


Phase 


Target 


J2000 


AOR key Date (JD) 


Class a 


(days) 


<J> <H> <K S > 


AK S Zero (JD) Obs. b Not 



NGC 362 V2 


01 


03 


21.85 


-70 


54 


20.1 


21740800 


2454268 


NGC 362 V16 


01 


03 


15.10 


-70 


50 


32.3 


21740800 


2454268 


NGC 5139 V42 


13 


26 


46.36 


-47 


29 


30.4 


21741056 


2454310 


NGC 5904 V84 


15 


18 


36.15 


+02 


04 


16.3 


21741312 


2454316 


NGC 5927 VI 


15 


28 


15.17 


-50 


38 


09.3 


21741568 


2454348 


NGC 5927 V3 


15 


28 


00.13 


-50 


10 


24.6 


21741824 


2454348 


Lynga 7 VI 


16 


11 


02.05 


-55 


19 


13.5 


21742080 


2454350 


NGC 6171 VI 


16 


32 


24.61 


-13 


12 


01.3 


21742336 


2454345 


NGC 6254 V2 


16 


57 


11.74 


-04 


03 


59.7 


21742592 


2454347 


NGC 6352 V5 


17 


25 


37.52 


-48 


22 


10.0 


21742848 


2454384 


NGC 6356 VI 


17 


23 


33.72 


-17 


49 


14.8 


21743104 


2454375 


NGC 6356 V3 


17 


23 


33.30 


-17 


48 


07.4 


21743104 


2454375 


NGC 6356 V4 


17 


23 


48.00 


-17 


48 


04.5 


21743104 


2454375 


NGC 6356 V5 


17 


23 


17.06 


-17 


46 


24.5 


21743104 


2454375 


NGC 6388 V3 


17 


36 


15.04 


-44 


43 


32.5 


21743360 


2454377 


NGC 6388 V4 


17 


35 


58.94 


-44 


43 


39.8 


21743360 


2454377 


Palomar 6 VI 


17 


43 


49.48 


-26 


15 


27.9 


21743616 


2454374 


Terzan 5 V2 


17 


47 


59.46 


-24 


47 


17.6 


21743872 


2454374 


Terzan 5 V5 


17 


48 


03.40 


-24 


46 


42.0 


21744128 


2454374 


Terzan 5 V6 


17 


48 


09.25 


-24 


47 


06.3 


21744128 


2454374 


Terzan 5 V7 


17 


47 


54.33 


-24 


49 


54.6 


21744128 


2454374 


Terzan 5 V8 


17 


48 


07.18 


-24 


46 


26.6 


21744128 


2454374 


Terzan 5 V9 


17 


18 


11.86 


-24 


50 


17.1 


21743872 


2454374 


NGC 6441 VI 


17 


50 


17.09 


-37 


03 


49.7 


21744384 


2454577 


NGC 6441 V2 


17 


50 


16.16 


-37 


02 


40.5 


21744384 


2454577 


NGC 6553 V4 


18 


09 


18.84 


-25 


54 


35.8 


21744640 


2454373 


IC 1276 VI 


18 


10 


51.55 


-07 


10 


54.5 


21742080 


2454376 


IC 1276 V3 


18 


10 


50.79 


-07 


13 


49.1 


21742080 


2454376 


Terzan 12 VI 


18 


12 


14.18 


-22 


43 


58.9 


21745152 


2454373 


NGC 6626 V17 


18 


24 


35.84 


-24 


53 


15.8 


21745408 


2454375 


NGC 6637 V4 


18 


31 


21.88 


-32 


22 


27.7 


21745664 


2454582 


NGC 6637 V5 


18 


31 


23.44 


-32 


20 


49.5 


21745664 


2454582 


NGC 6712 V2 


18 


53 


08.78 


-08 


41 


56.6 


21745920 


2454385 


NGC 6712 V7 


18 


52 


55.38 


-08 


42 


32.5 


21745920 


2454385 


NGC 6760 V3 


19 


11 


14.31 


+01 


01 


46.6 


21746176 


2454381 


NGC 6760 V4 
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a SR = semi-regular; CW = Cepheid of the W Vir type (i.e. Pop. II). 

b Phase during the IRS observation for P > 50 days. 

c Photometry from 2MASS, variability class and period from |Clcmcnt ctlil] (|2001l ). 



is particularly useful in the crowded regions typical for 
globular clusters. 

We used the optimal extraction for all of our data ex- 
cept for the LL portion of three targets: NGC 362 V16, 
Terzan 5 V5, and NGC 6441 V2. In these cases, the 
optimal algorithm was unable to separate our intended 
target from adjacent sources, and for these, we used the 
tapered-column extraction. 

The calibration of SL spectra is based on observations 
of HR 6348 (K0 III). The calibration of LL used HR 
6348 along with the late K giants HD 166780 and HD 
173511. We extracted the spectra of these stars using 
the optimal and tapered-column algorithms to calibrate 
the two extraction methods independently. 

The final step in the data reduction corrects for dis- 
continuities between the orders and trims untrustworthy 
data from the ends of each order. The correction for 
discontinuities applies scalar multiplicative adjustments 
to each order, shifting them upwards to the presumably 
best-centered segment. The exception is the three tar- 
gets blended with other sources in LL. For these, we nor- 



malized to SL1 to minimize the impact on our measured 
bolometric magnitudes of the contamination from the ad- 
ditional sources in LL. 

Figures [2] through [8] present the resulting spectra for 
the targets in Table [2 organized by the metallicity of 
the clusters of which they are members (and ordered by 
the right ascension of the cluster). Figure [7] includes 
the four Cepheids, three of which were observed with SL 
only, and Figure [S] includes those sources which are not 
oxygen- rich (Lynga 7 VI), are members of clusters with 
highly uncertain metallicities (Palomar 6 VI), or whose 
membership or evolutionary status is uncertain (§ 13 .T|) . 

2.4. Contemporaneous Photometry 

Table [3] presents photometry of 31 of our 39 targets 
obtained from the 2.3-m telescope at Siding Spring Ob- 
servatory (SSO) during the period when the Spitzer 
spectra were obtained. All observations were made us- 
ing the Cryogenic Arr ay Spectrometer/Imager (CASPIR; 
McGregor et al.lll994l) with the following filters: J (effec- 
tive wavelength 1.24 /im), H (1.68 /xm), K (2.22 /jm), and 
narrow-band L (3.59 /im). Calibrations were based on 
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Fig. 1. — Light curves in the K band for the sources in our sample observed with SIRIUS at the IRSF 1.4-m reflector. Table [2] presents 
the fitted periods, mean magnitudes, K-band amplitude, and zero-phase dates. 



observ ations of standard stars from the lists of lMcGregorl 
(1994), and the data were processed using standard tools 
available with IRAF (the Image Reduction and Analysis 
Facility) . 

3. MEMBERSHIP 

The principal objective of this study is to investigate 
how the quantity and composition of the dust produced 
by AGB stars depend on its initial metallicity. We as- 
sume that the metallicity of a star is simply the metal- 
licity of its cluster, making it important to identify non- 
members which are in the foreground or the background 
of a cluster. Our primary tools for testing membership 



are whether the relation between period and absolute K 
magnitude and the bolometric magnitude of the stars 
produces a distance consistent with membership, and if 
that distance results in a bolometric magnitude consis- 
tent with the masses of stars as old as those in our sam- 
ple. 

3.1. The Period-K Relation 

We estimated the distance to each v ariable by apply- 
ing th e period-K magnitude relation of iWhitelock et all 
(2008 ). In gen e ral, M K = p (log P - 2.38) + 5. 
IWhitelock et all ((2008) found that p = -3.51 in Mira 
variables in the LMC and the Galaxy, but 6 var- 
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Fig. 2. — IRS spectra of the five sources in our sample with 
metallicities ([Fe/H]) between 0.0 and —0.25. The error bars are 
generally smaller than the width of the plotted spectra. 

ied from -7.15 in the LMC to -7.25 in the Galaxy. 
iWhitelock et al.l ()20Q8f ) assumed that the distance mod- 
ulus to the LMC was 18.39. The value for 5 would be 
—7.25 for both the Galactic and LMC samples if the 
distance modulus to the LMC were 18.49 . whic h com - 
pares favorably to recent measurements. lAlvesI (|2004fi 
estimated the distance modulus to the LMC to be 
18.50±0.02, based on a revi e w of 14 recent determi- 
nations, and Keller fc Woo d (2006) found a value of 
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Fig. 3. — IRS spectra of the five sources in our sample with 
metallicities ([Fe/H]) between —0.25 and —0.49. 

18.54±0.02, based on an analysis of bump Cepheids. 
Thus, we assume that 5=— 7.25 for all sources in the 
following analysis. 

We dereddened the K magnitude using the estimates 
for E(B — V) given in Table Q] and the int e rstella r 
extinction measurements by iRieke &: Lebofskvl |l985), 
who find that R = A V /E{B - V) = 3.09 ± 0.03 and 
Ak/A v = 0.1123 We excluded the four Cepheid vari- 

16 We will also use the ratios Aj/A v = 0.282, A H /A V = 0.175, 
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Fig. 4. — IRS spectra of the nine sources in our sample with 
metallicities ([Fe/H]) between —0.49 and —0.65. 

ables in our sample from this analysis, which accounts 
for all of the periods less than 50 days. 

The period-K relation produces distance moduli for 
IC 1276 VI an d V3 of 13.51 and 13 .81, respective l y, while 
Harris I (|2003l) reports 13.66 and iBarbuv et all (|1998bO 
report 13.01. The latter is inconsistent with our esti- 
mates, and we adopt the former. Similarly, the period- 
K relation gives a distance m odulus for Te rzan 12 VI 
of 13.75, comp ared to 13.38 (jHarris I l200l and 12.66 
(jOrtolani et al.l fl998h . This source has a period of 458 
days, and we might expect circumstellar extinction at K 
of a few tenths of a magnitude. Thus we drop 12.66 as 
inconsistent with the period of the star and adopt 13.38 
as the distance modulus. 

Five of the six variables in Terzan 5 have an aver- 
age distance modulus of 14 .11 (±0. 13), compared to va l- 
ues of 15.0 5 (lHarris 1 [20031) . 13 87 (IValenti et alJl200l . 
and 13.70 (|Ortolani et all 120071) . Terzan 5 V9 is heav- 
ily reddened at K, giving an apparent distance modu- 
lus of 15.86, widely at variance with the others. Dif- 
ferential reddening in front of Terzan 5 would have to 
produce 15.6 additional magni tudes of visual extin ction 
and is an unlikely explanation. iFerraro et al.l (|2009l ) have 
identified a younger and more metal-rich population in 

and A^/Ay = 0.058 elsewhere in this paper. 
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Fig. 5. — IRS spectra of the six sources in our sample with metal- 
licities ([Fe/H]) from -0.65 to -0.80. 
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Fig. 6. — IRS spectra of the five AGB variables in our sample 
with metallicities ([Fe/H]) below -0.90. 

Terzan 5 (age^ 6 Gyr and [Fe/H]~+0.3). They present 
isochrones showing that the younger population should 
have roughly one extra magnitude of reddening in V— K, 
but their isochrones are based on models which do not 
extend to the tip of the AGB. We cannot rule out the pos- 
sibility that Terzan 5 V9 is a member of this younger pop- 
ulation, but it clearly differs from the remaining sources 
in the cluster and to be cautious, we will treat it as a 
non-member. The distance moduli of the remaining five 
variables are mutually consistent with each other and 
nearly two standard deviations from the closest published 
estimate, 13.87. For the remainder of this paper, we 
will adopt our own mean distance modulus for Terzan 5: 
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Fig. 7. — IRS spectra of the four Cepheid variables in our sample. 
All have metallicities ([Fe/H]) below —1.0. Only one of the four 
has LL data. 
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Fig. 8. — IRS spectra of the five sources which cannot be in cluded 
in the metallicity-dependent analysis, as explained in § I3.3I 

14.11. 

Figured] presents the period-K relation for our sample 
of globular cluster variables (excluding the Cepheids). 
The absolute K magnitudes are based on the distance 
moduli presented in Table [T] The labelled sources in 
the figure require special consideration, which in some 
cases has led to their exclusion as members. In Fig- 
ure [9j though, their absolute magnitude is determined 
as though they were members. The figure also shows the 
period-K relation for fundamental-mode pulsators and 
for overtone pulsators, assuming that their periods are 
2.3 times smaller. 

Two sources appear to be too bright for the period-K 
relation: NGC 6171 VI and Palomar 10 V2. While it 
is possible that they are overtone pulsators, their bolo- 
metric magnitudes are inconsistent with membership, as 
explained in the next section. Assuming that both are 
fundamental-mode pulsators, NGC 6171 VI has a dis- 
tance modulus of 12.15, compared to 13.89 for the clus- 
ter, putting it 3.3 kpc in front of the cluster and 1.1 kpc 
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Fig. 9. — The periods and absolute K magnitudes of the vari- 
ables in our globular cluster sample (excluding Cepheids). Abso- 
lute magnitudes are determined from the distance moduli for the 
clusters reported in Table [TJ The solid line is the nominal period-K 
relation defined by Whitelock et al. (2008). The dashed line rep- 
resents the shift in the relation expected for stars pulsating in the 
first overtone mode. Sou rces which do not follow the period-K 
relation are labelled (see § 13.41 1. 

above the Galactic plane. Similarly, Palomar 10 V2 has a 
distance modulus of 12.30 vs. 13.86 for the cluster, which 
means it is 3.0 kpc closer and about 0.14 kpc above the 
Galactic plane. 

The three sources below the period-K relation present 
more complex cases. We will examine them after deter- 
mining bolometric magnitudes (next section). 

Figure |9] also shows that all of the confirmed members 
are pulsating in the fundamental mode, even though sev- 
eral sources are classified as semi-regular variables, or 
possible semi- regular variables, in Table [5J These results 
are consistent if they are SRa variables. 

3.2. Bolometric Magnitudes 

Table |4] includes bolometric magnitudes for each star 
in our sample. Basically, we integrated the photome- 
try from SSO and the spectroscopic data from Spitzer 
and corrected for the distances in Table Q] (except for 
the two most likely non-members, NGC 6171 VI and 
Palomar 10 V2, where we used the distance moduli de- 
termined in the preceding section). For the eight stars 
not observed from SSO, we substituted the mean mag- 
nitudes in Table [5] Most of the SSO photometry was 
taken within one or two days of the Spitzer observations, 
but when this was not the case (and for the non-SSO 
photometry), we phase-corrected the data to match the 
Spitzer epoch, using the amplitudes from our analysis 
of the IRSF/SIRIUS monitoring data. The monitoring 
data did not include observations at L, and we assumed 
that AL~ AK, as confirmed by the mean amp litudes for 
oxygen-rich LPVs published by iSmithl (|2003l ). We also 
corrected the photometry for interstellar extinction, us- 
ing the E(B — V) excesses in Tab le [Hand the insterstellar 
reddening of iRieke fc Lebofskvl (|1985| ). To integrate the 
flux outside of the observed wavelengths, we extended a 
3600 K blackbody to the blue and a Rayleigh- Jeans tail 
to the red. Finally, to correct the bolometric magnitudes 
for pulsation, we assumed that the K-band amplitude ap- 
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TABLE 4 

Spectroscopic Properties 



Target 


M bol 




[7] 


-[15] 




DEC a 


lop 

(M 


; M 

/yr) 


b 


Corrected 
Fu/F 12 


IR Spec. 
Class c 


Notes 


IN GO oOz VZ 


—3.50 





.lb 


± 





.02 





.09 


± 


0.01 


—7.90 


± 


0. 


.06 




l.XN 






IN GO ooz Vlb 


—4.10 


0. 


.25 


± 


0. 


.02 


0. 


.02 


± 


0.00 


—8.54 


± 


0. 


11 




l.N 






INGG 51o9 V4z 


—3.88 





.41 


± 





.01 


0. 


07 


± 


0.00 


—8.02 


± 


0. 


.11 




l.N 






INGG oyU4 Vo4 


—3.05 












— 0. 


06 


± 


0.03 












1 TVT 
1.1N 


d, 


c 


IN GO 59z7 VI 


—4.13 


1 


.19 


± 


0. 


.02 


0. 


67 


± 


0.01 


—6.65 


± 


0. 


11 


1 00 1 n m 

1.22 ± U.U1 


2.SX4t 






in go oyz < v o 


— 4.64 


1 


.35 


± 





.00 


1. 


.48 


± 


0.00 


—6.27 


± 


0. 


03 


1 c:7 _i_ n n c; 

1.57 ± U.Uo 


2.SE8t 






Lynga 7 VI 


—4.99 












0. 


58 


± 


0.01 


—6.61 


± 


0. 


.01 




Z.Oil/ 






IN GO 01 f 1 VI 


—5.11 


1 


.53 


± 





.01 


1. 


.99 


± 


0.01 


—6.02 


± 


0. 


.07 


l.bo ± U.Uo 


2.ohio 


f 




IN GO ozo4 Vz 


— 1.99 












— 


02 


± 


0.01 












l.N 




c 


INGG oooz Vo 


—4.27 


1 


.05 


± 





.02 





.59 


± 


0.01 


—6.81 


± 


0. 


.01 


1 ie _i_ n no 
l.lo ± U.Uo 


2.o A4t 






IN GO oooo V 1 


— 4.54 


1 


.1/ 


± 





.01 


2. 


08 


± 


0.01 


—6.06 


± 


0. 


.03 


l.OZ ± U.Uo 


Z.bil/or 






■\t/~i ^ core i ro 
IN GO OOOO VO 


—4.06 





.58 


± 





.02 





17 


± 


0.01 


—7.66 


± 


0. 


02 


r\ no 1 r\ r\ A 

0.83 ± 0.04 


2.SE1 






IN GO boob V4 


—4.42 





.61 


± 





.02 


0. 


.14 


± 


0.01 


—7.57 


± 


0. 


.03 


1 no 1 n n cr 

l.Uo ± U.Uo 


2.SE2 t: 






IN GO oooo Vo 


— 4.61 





p. ~ 
.U t 


± 





.02 


0. 


25 


± 
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—7.40 


± 


0. 
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l.Oi ± U.U4 
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—3.83 





.13 
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.02 
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l.N 






IN GO oo8» V4 
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68 
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± 
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1 to 1 n r\ r ? 

1.73 ± 0.07 


2.SE8 






Falomar o VI 


—3.96 


2 


.02 
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.01 


1 


27 


± 


0.01 


—5.11 


± 


0. 


.43 




3.SBxf 


g 




Terzan 5 V2 


—4.41 





.62 


± 





.01 


0. 


11 


± 


0.01 


—7.70 


± 


0. 


17 


0.79 ± 0.07 


2.SE1 t: 






lerzan 5 Vo 


—5.09 


1 


.38 


± 





.02 


0. 


.91 


± 


0.02 


—6.39 


± 


0. 


22 


1 00 1 n no 

1.23 ± 0.02 


2.SE4 






lerzan 5 Vo 


—4.72 


1. 


.60 


± 





.01 


1. 


11 


± 


0.01 


—6.06 


± 


0. 


30 


1 1 T 1 n m 

1.17 ± U.U1 


2.SX4t 






1 .. ...... ... C ^ T*7 

lerzan o V / 


—4.98 


1 


.61 


± 





.01 


1. 


71 


± 


0.00 


—6.00 


± 


0. 


23 


l.oz ± U.Uo 


Z.oH/O 






lerzan o Vo 


—5.04 





.78 


± 





.01 


0. 


15 


± 


0.01 


—7.46 


± 


0. 


21 


a 01 _i_ nil 
U.zl zb U.ll 


Z. oil/1 t: 






lerzan o Vy 


—3.72 


1 


.46 


± 





.03 


0. 


.88 


± 


0.01 


—6.33 


± 


0. 


33 


n nt: j_ n no 

U.yo ± U.Uo 


O CVO+ 

2.o 1 2t 


g 




In LrO 0441 V 1 


— 4.26 





.71 


± 





.01 


0. 


20 


± 


0.01 


— 7.44 


± 


0. 


03 


n A _L n n a 
U.o4 ± U.U4 


O QTPI 4- 






JNtjO 6441 V2 


—3.98 


0. 


.52 


± 


0. 


.03 


0. 


.12 


± 


0.02 


—7.76 


± 


0. 


.01 


n nn I n t a 

0.99 ± 0.14 


2.SE2 






1N(jO 6553 V4 


—4.60 





.90 


± 





.01 





.32 


± 


0.01 


—8.60 


± 


0. 


57 


n n A 1 n nn 

0.94 ± 0.02 


2.SYlt 






lO 1-2 iu VI 


—4.60 





.71 


± 





.01 





13 


± 


0.01 


—7.57 


± 


0. 


21 


0.b8 ± 0.08 


O CT7 1 1 +. 

Z. oil/1 t: 






IC 1276 V3 


-4.72 


1 


.02 


± 





.02 


0. 


.52 


± 


0.01 


-6.88 


± 


0. 


.01 


1.22 ± 0.01 


2.SE4 






Terzan 12 VI 


-4.72 


1 


.48 


± 





.02 


1 


06 


± 


0.01 


-6.26 


± 


0. 


27 


1.13 ± 0.02 


2.SY3 






NGC 6626 V17 


-3.42 





.27 


± 





.02 


-0. 


.08 


± 


0.01 


-8.60 


± 


0. 


.57 




l.N 


d 




NGC 6637 V4 


-4.44 


1 


.06 


± 





.01 


1. 


.12 


± 


0.02 


-6.61 


± 


0. 


27 


1.70 ± 0.08 


2.SE8 






NGC 6637 V5 


-4.44 





.56 


± 





.01 


0. 


15 


± 


0.01 


-7.66 


± 


0. 


03 


1.17 ± 0.04 


2.SE4 






NGC 6712 V2 


-3.91 





.56 


± 





.02 


0. 


27 


± 


0.01 


-7.48 


± 


0. 


.29 


1.62 ± 0.08 


2.SE8 






NGC 6712 V7 


-4.34 





.52 


± 





.01 





21 


± 


0.01 


-7.59 


± 


0. 


23 


1.51 ± 0.05 


2.SE7 






NGC 6760 V3 


-4.71 





.78 


± 





.01 


0. 


.51 


± 


0.02 


-7.09 


± 


0. 


.28 


1.72 ± 0.07 


2.SE8 






NGC 6760 V4 


-4.34 





.96 


± 





.01 


0. 


.60 


± 


0.01 


-6.88 


± 


0. 


13 


1.33 ± 0.02 


2.SE5 






NGC 6779 V6 


-2.96 












-0. 


06 


± 


0.01 












l.N 


d, 


c 


Palomar 10 V2 


-4.96 


1 


.21 


± 





.02 





68 


± 


0.00 


-6.63 


± 


0. 


.13 


1.14 ± 0.02 


2.SE3t 


h 




NGC 6838 VI 


-4.00 





.77 


± 





.01 


0. 


.36 


± 


0.01 


-7.21 


± 


0. 


.15 


1.42 ± 0.04 


2.SE6 t: 







a Dust emission contrast; see §4.1 for an explanation and a discussion of systematic errors. 

b Subtract 2.30 to convert to login of the dust-production rate (see §4.2). 

c These classifications are defined in §4.1. 

d Cepheid variable. 

c SL data only. 

f Assumed m - M = 12.15. 

g Membership uncertain. 

h Assumed m - M = 12.30. 



proximated the bolometric amplitude and used the phase 
to correct to the mean. 

Figure [TO] plots the bolometric magnitudes of our 
sample against their pulsation periods, and it in- 
cludes (simplified) evolu tionary tracks taken from 
IVassiliadis fe Wood! (Il99l Fig. 20). We will return to 
this diagram in §[6] Here, the diagram can help us assess 
cluster membership. Except for two sources in the lower 
right (Terzan 5 V9 and Palomar 6 VI), the entire sample 
appears to follow a single sequence from the lower left to 
the upper right. Terzan 12 VI is the orange diamond 
at the lower edge of the sequence with a period of 458 
days. It is close enough to the sequence defined by the 
remaining sources that it validates our choice of distance 
modulus. The next section examines the three sources 
that appear below the period-K relation (Figure ^ in 
turn. 



3.3. Special Cases 

The bolometric magnitudes of NGC 6171 VI and Palo- 
mar 10 V2 in Table [4] are based on the presumption 
that they are foreground objects (§ 13. ip . If they were 
actually overtone pulsators and cluster members, their 
bolometric magnitudes would be —6.85 and —6.52, re- 
spectively. Accord i ng to the evolutionary models of 
IVassiliadis fc Wood! (|l99l . these magnitudes would cor- 
respond to initial masses of ^5 M Q , much too massive 
for the age of either cluster. 

Terzan 5 V9 is about 1.7 magnitudes too faint at K if 
it were at the distance of the cluster. Its J— K and H— K 
colors, corrected for interstellar extinction, are consistent 
with 1.2-1.3 magnitudes of circumstellar extinction at K, 
and its K— L color suggests A/f~1.6. These color-based 
estimates assume that the circumstell ar dust extinction 
follows the interstellar relationships of lRieke fc Lebofskyl 
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Fig. 10. — The bolometric magnitudes and pulsation periods of 
the globular sample, compar ed to models of AGB sta rs at 0.945, 
1.0, and 1.5 M Q (black lines; IVassiliadis fc W ood 1993;). The sym- 
bols are color-coded to run from the most metal-rich globular clus- 
ter (red) to the most metal-poor (blue), with special cases in brown. 



(|1985ft . and they are roughly consistent with the observed 
extinction. However, this source sits in the lower right in 
Figure [10] with P=464 days. If it really is a member of 
Terzan 5, its bolometric magnitude is inconsistent with 
its pulsation period. Perhaps it is a background object, 
or perhaps it is a binary or interacting system. Whatever 
it is, it is not a normal AGB star, and we consequently 
do not include it when analyzing the other AGB stars in 
Terzan 5. 

Palomar 6 VI is too faint by 2.6 magnitudes at K, 
but its J— K, H— K, and K— L colors correspond to cir- 
cumstellar extinction at K of 3.5-4.9 magnitudes. These 
estimates would put the source ~l-2 magnitudes of dis- 
tance modulus in front of Palomar 6. In Figure [101 Palo- 
mar 6 VI is the source below the others in the lower 
right with P=566 days. If it were in the foreground of 
Palomar 6, correcting for its distance would push it even 
further down in Figure [TU] Here again, the star cannot 
be a normal AGB star and a cluster member. As already 
noted f ^2.ip . even if it were, the uncertain metallicity of 
the cluster prevents us from assigning it to one of our 
metallicity bins. 

The spectrum of Lynga 7 VI (Figure [8J is unambigu- 
ously that of a carbon star, with strong dust emission 
from SiC at ~11.5 fxm and MgS dust at ~26 fim, as well 
as absorption from acetylene gas at 7.5 and 13.7 fxm. Its 
carbon-rich status may be consistent with the circumstel- 
lar extinction at K apparent in Figure |9l but it would be 
helpful to independently verify its di stance. 

We have two means of doing so. iSloan etlrtl J2008) 
calibrated a relationship of Mk vs. J— K color for carbon 
stars in the SMC, using 2MASS photometry Correcting 
the mean magnitudes in Table [2] for interstellar extinc- 
tion give s J— K = 3.71 , whic h leads to a distance modulus 
of 14.71. ISloan et al.l ()2008j ) also found and calibrated a 
color-magnitude relation for the narrow 6.4- and 9.3-/Ltm 
filters defined to analyze earlier samples of Magellanic 
carbon stars. Lynga 7 VI has a [6.4] — [9.3] color of 0.55, 
which implies M 9 . 3 =-11.59. Since [9.3]=2.67, the dis- 
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Fig. 11. — A comparison of the bolometric magnitudes of our 
globular sample o f AGB stars and the the Magellanic sample of 
ISloan et al.l ( 2008). The histogram of the globular sample excludes 
the four Cepheid variables (all with bolometric magnitudes between 
— 1.6 and —2.7) and the four sources whose cluster membership is 
in doubt. 

tance modulus is 14.26. The average from these two 
methods is 14.20 ± 0.08, which compares favorably to 
the nominal distance modulus of Lynga 7, 14.31F1. 

In Figure [TQ1 Lynga 7 VI is the upper right datum. Its 
position is consistent with the evolutionary sequence de- 
fined by the larger sample, adding some confidence that it 
is a cluster member. Nonetheless, its carbon-rich nature 
prevents any comparison with the oxygen-rich sample, 
and it remains excluded from most of the analysis in this 
paper. 

3.4. Membership Summary 

To summarize, we exclude NGC 6171 VI and Palo- 
mar 10 V2 from the sample because if they were cluster 
members, they would be too bright for old AGB stars. 
We exclude Terzan 5 V9 and Palomar 6 VI from fur- 
ther consideration because if they are cluster members, 
they are too faint for their pulsation periods compared 
to the rest of the sample. Finally, while Lynga 7 VI is a 
cluster member, it is a carbon star and must be treated 
separately from the rest of the sample. 

Figure [TT] compares the bolometric magnitudes of the 
31 LPVs in our sample which are confirmed as members 
of globular clusters and normal AGB stars to the oxygen - 
rich Magellanic sample considered by Slo an et all ((2008) . 
The Magellanic sample spans a narrower range of metal- 
licity (—0.6 < [Fe/H] < —0.3). It is readily apparent 
that the samples represent very different sources. The 
Magellanic sample consists primarily of supergiants and 
brighter AGB stars, while the globular sample contains 
only of fainter and lower-mass AGB stars. This differ- 
ence arises primarily from the greater distance to the 
Magellanic Clouds and the resulting bias toward more 
luminous sources. 

4. SPECTRAL ANALYSIS 
4.1. Infrared Spectral Classification 

17 Using an u pdated calibration o f Mj( vs. J— K for the more 
metal-rich LMC (Lagadcc et al. 2010) gives m — M = 14.71 and a 
mean m — M = 14.49 ± 0.32, which is still close to the nominal 
distance. 
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Fig. 12. — The continuum-subtracted flux ratios of the globular 
sample compared to the silicate d ust s equence (dashed line). The 
symbols are color-coded as in Fig. 1101 The circled brown cross in 
the lower left is the carbon star Lynga 7 VI. The four points to the 
right of the silicate dust sequence (with 1.3 < F\o/Fvz < 1-8) are, 
from top to bottom, NGC 6712 V7, NGC 6838 VI, NGC 6356 V5, 
and NGC 6760 V4, and their deviation from the silicate dust se- 
quence may result from as-yet-unidentified impurities in the dust. 



The classificatio n of the spectra follo ws the Hanscom 
system, defined bv lKraemer et all (|2002fl for spectra from 
the Short- Wavelength Spectrometer (SWS) aboard the 
Infrared Space Observatory (ISO), which is based par- 
tially on the classificatio n of spectra of o xygen -rich AGB 
variables developed by iSloan fc Pried (|1995l) for data 
from the Low-Resolution S pectrograph o n the Infrared 
Astronomical Satellite. Sloa n et al.l (|2008f) explained how 
the method was modified for IRS data from Spitzer. In 
the Hanscom system, spectra are divided into groups, 
based on their overall color, and one- or two-letter des- 
ignations are added to describe the dominant spectral 
features. All of our spectra fall into Group 1 (for blue 
spectra dominated by stellar continua and showing no ob- 
vious dust), Group 2 (for stars with dust), and Group 3 
(for spectra dominated by warm dust emission) . Most of 
our spectra can be classified as either naked stars ( "l.N" ) 
or stars showing silicate emission ("2.SE."). 

The classification depends primarily on two quantities. 
The first quantity is Dust Emission Contrast (DEC), de- 



fined as the ratio of the dust excess to the stellar contin- 
uum, integrated from 7.67 to 14.03 /im. For the stellar 
continuum, we assume a 36 00 K Planck function , fit to 
the spectrum at 6.8-7.4 /im. ISloan fe Price! (fl99l 1998) 
assumed an Engelkc function with 15% SiO absorption at 
8 /xm. Switching to this continuum would systematically 
shift our DEC measurements upward by ^0.04. Table 0] 
does not include this systematic error in the uncertain- 
ties, although it should be kept in mind that the actual 
stellar continuum could vary from one source to the next. 
In this sample and with these assumptions, a DEC of 0.10 
separates those stars which we visually identify as naked 
from those with apparent dust excesses. 

The second quantity, determined for the spectra with 
an oxygen-rich dust excess, is the ratio of the excess emis- 
sion at 11 and 12 /im (F n/J 7 ^)- To measur e this fl ux ra- 
tio, we follow the method of ISloan fc Pried (]1995f ). mea- 
suring the excess at 10, 11, and 12 /im and plotting the 
flux ratio F w /Fi2 as a function of F w /Fu. Figure [T2l 
shows that all oxygen-rich sources fall on or close to the 
silicate dust sequence, which was defined as a power law: 
F w /F 12 = 1.32 (Fw/Fu) 1 - 77 . In the bottom left, the 
point furthest from the power law is the carbon star 
Lynga 7 VI. For the remaining sources, the flux ratio 
F\i/Fi2 is determined by finding the point on the sil- 
icate dust sequence closest to the point defined by the 
measured ratios fio/^ii and -Fio/-Fi2- The corrected ra- 
tio defines the silicate emission (SE) index, which runs 
from 1 for Fu/F 12 < 0.85 to 8 for Fu/F 12 > 1.55. Four 
sources have positions shifted to the right of the sili- 
cate dust sequence in the region where 1.3 < F w /Fi 2 
< 1.8; from top to bottom, they are NGC 6712 V7, 
NGC 6838 VI, NGC 6356 V5, and NGC 6760 V4. Their 
spectra have typical silicate emission features, and we 
suspect that the shift arises from a dust component in 
addition to the usual mixture of amorphous alumina and 
amorphous silicates. 

The correc ted flux ratio FjjJF iz quantifies the dust 
composition (|Egan fc Sloanll2~001[ ). Amorphous alumina 
dominates the spectra with low flux ratios (or SE indices 
1-3), while amorphous silicates dominate the highest flux 
ratios (SE6-8). For the remainder of the paper, we will 
use the quanitity FwfFyi to distinguish the spectra by 
the composition of their dust. 

Lynga 7 VI is classified in the Hanscom system as 
"2.CE" (carbon-rich emission) . Of the remainder, 10 are 
naked ("l.N"), and 28 have oxygen-rich dust. Most are 
classified in the sequence from to "2.SE1" to "2.SE8", 
but seven spectra require special attention, as explained 
next. 

Strong silicate self-absorption can push spectra down 
the silicate dust sequence, as has happened for Palo- 
mar 6 VI. This source mimics a 2.SE1 spectrum in 
Figure but its spectrum (Figure [5]) clearly shows 
self-absorption at 10 /im. That, and the fact that the 
spectrum peaks past 15 /im, result in a classification of 
"3.SB". 

Six of the spectra are distinctly different from the rest 
and are examined more closely in § 14.51 and 14.61 below. 
Three clearly show the presence of crystalline silicates in 
the 10-/im feature; they are classified as "SX" instead of 
"SE". Three others show an unusual spectral emission 
feature peaking between 11 and 12 /im. To distinguish 
them from the remainder, we will give them the new (and 
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possibly temporary) classification of "SY" . Both the SX 
and SY sources are indexed identically to the SE sources, 
although the interpretation of that index may differ. 

The classifications in Table 0] include the suffixes 't', 
'£', and 'x', indicating the presence of a 13-/im feature, 
a 14-/im feature, and longer-wavelength features from 
crystalline silicates, respectively. These emission features 
and the criteria for their classification are treated in more 
detail below (§ S3] and |4jj) . 

4.2. Mass-loss Rates 

The primary objective of this project is to quantify 
how the rate of dust production by evolved oxygen-rich 
stars depends on metallicity. The dust-production rate 
is identical to the dust mass-loss rate, which is the ratio 
of the overall mass-loss rate divided by the gas-to-dust 
ratio. We use two methods to quantify the amount of 
dust in the circumstellar shell: t he DEC described a bove 
and the [7] -[15] color defined by ISloan et all (|2008D and 
described below. To relate these to the dust mass-loss 
rate, we use the set of radiative transfer models fit- 
ted to 86 evolved s tars in the Magellanic Clouds by 
iGroenewegen et aTl (|2009l ). Their models determine a 
dust mass from the opacity needed to fit the IRS spec- 
troscopy and overall spectral energy distribution, assume 
a constant outflow velocity (10 km s _1 ) to determine the 
dust mass-loss rate, and assume a gas-to-dust ratio of 
200 to determine the overall mass-loss rate. 

The [7] — [15] color integrates the total flux density 
(from star and dust) in the wavelength i nterva ls 6.8-7.4 
and 14.4-15.0 /im. IGroenewegen et ail (|2009f) recently 
showed that the [7] — [15] color tracks the mass- loss rate: 

logM(M /yr) = 1.759([7] - [15]) - 8.664. (1) 

They also found that the DEC correlates with mass-loss 
rate as well: 

log M(M /yr) = 1.392DEC - 6.484, (2) 

for mass-loss rates less than 10~ 55 M Q /yr. At higher 
rates, the shell becomes optically thick enough to drive 
the 10 /im silicate emission feature into self-absorption, 
and as a result, the DEC breaks down as a useful measure 
of dust content. 

For all sources except the carbon star Lynga 7 VI 
(§ 15.11) and three Cepheids not observed at 15 /im, we 
estimated mass-loss rates using Equations (1) and (2). 
We began with [7] — [15] color. For those sources where 
[7] — [15] > 1.80, indicating a mass- loss rate above 10 -55 
MQ/yr, we used only Equation (1). For those sources 
below this limit, we also estimated a mass- loss rate from 
the DEC and averaged the two. We assumed a minimum 
DEC of 0.0 for the purposes of estimating mass- loss rates. 

Table |4] presents the results as overall mass-loss rates, 
in order to facilitate comparison to other published mass- 
loss rates. To convert overall mass-loss rate to dust- 
production rate, one needs to divide by the assumed gas- 
to-dust ratio of 200, or in log space, to subtract 2.30. For 
the remainder of the paper, we will focus on the dust- 
production rate, which avoids assuming a gas-to-dust ra- 
tio. 

4.3. Dust Production Dependencies 



Generally, the mass-loss rate and dust production rate 
will increase as a star evolves up the AGB, due to the 
reduced surface gravity and higher luminosity of the cen- 
tral star. This increase will dominate the subtler effect of 
metallicity. To examine the role of metallicity, one would 
ideally plot dust production vs. luminosity for samples of 
different metallicities and compare them. Past compar- 
isons between evolved stars in the Galaxy and the Mag- 
ellanic Clouds have been hampered by the poorly con- 
strained distances to Galactic sources, and consequently, 
pulsation period has been used as a proxy for luminosity. 
Here, we can overcome this difficulty, since the distances 
to the globular clusters in our sample are known and we 
can determine bolometric magnitudes directly. 

Figure [13] plots the derived mass-loss rates as a func- 
tion of pulsation period (Panel a), bolometric magnitude 
(Panel b), and the corrected flux ratio Fn/Fi2, which 
quantitifies the dust composition (Panel c). In general, 
the rate of dust production clearly increases as the stars 
grow brighter, their pulsation periods increase, or the 
dust grows more silicate ri c h. Pa nel c is analogous to 
Figure 4 by iSloan fc Price) (|1995f ). which showed that 
alumina-rich dust (SE1-3) is seen only in low-contrast 
shells, while silicate-dominated dust (SE6-8) can exhibit 
a wide range of dust emission contrasts. The globular 
sample follows the same trend. 

To examine how the metallicity might influence the 
tendency of the dust production to increase with increas- 
ing period, luminosity, or silicate/alumina dust ratio, 
we have determined a mean slope for each metallicity- 
dcfincd subsample, including all the spectra depicted in 
Figures [5] through assuming that they have same y- 
intercept (which we arbitrarily defined). These slopes 
appear in Figure [13] as dotted lines. Figure [14] plots the 
slopes of these fitted lines as a function of metallicity, 
showing how the metallicity modifies the more obvious 
relations of dust production with period, luminosity, and 
dust content. The uncertainties in Figure [TJ] are the for- 
mal uncertainties in the mean (the standar d deviation 
divid ed by the square root of the sample size; iBevington] 
[l969h . 

Panel a in Figures [13] and [14] shows the effect of metal- 
licity when segregating the samples by period, or more 
to the point, the lack of an effect. While the data may 
appear to show a negative relation, a horizontal line in 
Panel a of Figure [14] can touch all of the error bars. This 
lack of a dependence on met allicity is a little suprising 
given that ISloan et "all (|2008T l found a metallicity depen- 
dence in dust production when segregating by period and 
comparing oxygen-rich evolved stars in the Galaxy and 
the Magellanic Clouds. However, the globular sample 
spans a wider range of metallicities, and the dependence 
of pulsation period on metallicity may be obscuring the 
dependence of dust production on metallicity. One can 
see in Figure Q2] that no star in the most metal-poor 
bin has a period greater than 200 days, while no star in 
the most metal-rich bin has a period less than 20 days. 
The in termediate bins have intermediate periods. iWoodl 
(1990) noted that lowering the metallicity of an AGB 
star raises its effective temperature and thus reduces its 
radius, which also reduces its pulsation period. 

Panel b in Figures IT31 and [T4l shows that segregating by 
bolometric magnitude reveals a slight dependence of dust 
production rates on metallicity, although it is blurred 
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Fig. 13. — The log ol the dust-production rate (DPR) for each of the stars in our sample as a function of pulsation period (Panel a), 
bolometric magnitude (Panel b), and corrected flux ratio F11/F12 (Panel c). The overall mass-loss rate is the product of the DPR and 
the gas-to-dust ratio, which in Table [4] is assumed to be 200 ( a di fference of 2.30 in log space). The corrected flux ratio in Panel c locates 
the position of a source along the silicate dust sequence (Fig. I12H and measures the relative contributions from amorphous alumina (low 
F11/F12) and amorphous silicates (high F11/F12). The dotted lines give the mean slope of each metallicity group, assuming a common 
y-intercept, as explained in the text. 

somewhat by scatter in the sample. 

The impact of metallicity is clearest when segregat- 
ing the samples by dust content as measured by F\\jF\i 
(Panel c in the two figures). Whether the dust is 
alumina-rich or silicate-rich, the dust-production rate in- 
creases as the metallicity increases. 

Combining the evidence from the three methods of seg- 
regating the globular sample, we conclude that metallic- 
ity is influencing the rate of dust production by AGB 
stars. The effect may be stronger than what we ob- 
serve. Implicit in our calibration of dust-production 
rate from [7] — [15] and DEC is the assumption by 
IGroenewegen et all (|2009l ) that the outflow velocity is 10 
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Fig. 14. — The slopes from Fig. [13] plotted as a function of metal- 
licity, with vertical error bars based on the uncertainty in the mean 
(standard deviation divided by the square root of the sample size. 
The slopes are identical for overall mass-loss rate (MLR) and dust- 
production rate (DPR). These plots show more clearly how the 
mass-loss rate depends on bolometric magnitude and position along 
the silicate dust sequence (F11/F12), but not on pulsation period 
(as indicated by the horizon tal dott ed li ne in Panel a). Symbol 
shapes are as defined for Fig. 1101 1121 and 1131 



Galactic Halo by [Lagadec et al.l (|2010D reveal a possible 
dependence of the outflow velocity on metallicity. While 
this result requires confirmation, it is reasonable to con- 
sider the possibility that the outflow velocity increases 
with metallicity in our sample as well. In this case, we 
would have to revise our dust-production rates upward 
for the higher metallicities, and the trends in Fig.fTSIand 
1141 would be more apparent. 

On the other hand, if the gas-to-dust ratio increases 
at lower metallicity, this would decrease the dependency 
for overall mass-loss rate. The likelihood depends on how 
important the role of dust is in the mass- loss process from 
the AGB. 

4.4. Amorphous alumina 

Figure and Figure Q2] (Panel c) show that all of 
the metallicity bins contain spectra classified as SE1- 
3, whi ch arise from shells dominated by amorphous alu- 
mina (|Egan fc S loan 200l|). We conclude that AGB stars 
can produce alumina-rich dust, regardless of their initial 
metallicity. 

Sloa n et "ail (|2008T ) found little incidence of alumina- 
rich dust in their Magellanic samples, which they at- 
tributed to a possible underabundance of aluminum in 
more primitive stars, but our detection of alumina-rich 
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Fig. 15. — The three spectra classified as "SX" due to the likely 
presence of crystalline silicate emission at 10 and 11 fim. These 
spectra are plotted after a stellar continuum has been fitted from 
6.8 to 7.4 fim and removed. All three show features at each of the 
wavelengths marked by vertical lines: 9.7, 11.3, 13, and 20 fim. 
Terzan 5 V6, with the "9.7"-/mi feature shifted to 9.3 /im and a 
clear contribution from silicates at 18 /an feature, may be affected 
by some self-absorption in the 10 fim silicate feature. 



TABLE 5 

Narrow dust emission feature strengths 

Feature strength/dust strength (%) a 
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a Ratio of feature strength to total dust emission 
from 5 to 35 fim; see §4.4. 



dust shell s at even lower metal licities repudiates that ar- 
gument. iSloan fc Pried (|1998| ) found that Galactic su- 
pergiants rarely produced alumina-rich dust shells. Fig- 
ure [TT1 shows tha t the Magellanic sample studied by 
Sloa net all £2008) is more luminous, and it follows that 
it generally contains more massive stars. The lack of 
low-mass stars in the Magellanic samples may explain 
the missing alumina-rich dust. 

4.5. Warm crystalline dust 

Figure [15] shows the three spectra classified as "SX" , 
based on the splitting of the 10 /im-emission feature from 
amorphous sili cate grains i nto tw o components at 9.7 
and 11.3 /im. ISloan et al.l (pOOfih found one source in 
the LMC which showed similar structure at 10 fim, HV 
2310, and they showed that an increased fraction of crys- 
talline silicate grains could explain the spectrum. They 
suggested that grains might form with more crystalline 
structure in l ow-density dust-fo rmation zones, as de- 
scribed below. ISloan et al.1 ([2008) added a second source 
in the LMC: HV 12667. 

All three globular SX spectra also show a strong 13- 
fim feature, as well as an additional component at 20 fim. 
These features do not appear in the spectra of HV 2310 
or HV 12667 in the LMC. However, at least 24 other 
spectra in the globular sample show at least one of these 
features. Table [5] presents the strengths of the features in 
the spectra where at least one was detected, expressed as 
a percentage of the total dust emission from 5 to 35 fim. 
We measured these strengths by fitting a line segment 
over the wavelength ranges given in Table [6] and integrat- 
ing in between. We repeated the process using a spline 



TABLE 6 

Wavelength intervals for extracting 
narrow dust emission features 



Feature (fim) 


Continuum intervals (fim) 


13 


12.20-12.35 13.50-13.65 


20 


18.40-18.80 21.10-21.60 



to estimate the (dust) continuum, integrating over the 
same wavelength range. Where both methods produced 
a continuum that followed the actual data to cither side 
of the feature, we averaged the result. The uncertainties 
in Table [5] are the larger of the propagated error in the 
t wo extraction s or th e standard deviation between them. 

ISloan et "all (|2003l ) found a correlation between the 
strength of the 13- and 20-fj.m features, quoting a Pear- 
son correlation coefficient of 0.82. We have calculated 
the same coefficient for all of the sources in the present 
sample to be 0.62, notably but understandably less given 
the lower S/N and spectral resolution of the current set 
of spectra. 

Following [Sloan et al.l (|2003l) . we assign a suffix of "t" 
to the classification of all sources with a 13 fim feature 
stronger than 0.1% of the total dust emission, provided 
that the S/N of the detection is 3.0 or more. Where the 
feature exceeds 0.1% of the total dust, but the S/N is 
1.0-3.0, we classify the spectrum as "t:". Eight spectra 
have clear 13-/xm features, and six more have probable 
13-fim features, which accounts for 14 out of 30 spectra 
showing oxygen-rich dust emission, or 47%. Nearly all 
of the variables in the sample are classified as Miras, 
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and Figure [S] reinforces this point. Slo an et all (|1996ft 
estimated that only ~20% of Galactic Miras with oxygen- 
rich dust showed 13-/tm features; the percentage here is 
more than twice as high. This difference might be due 
to the lower masses in the globular sample. 

The specific grain which produces the 13-/im fea- 
ture remains an unsettled issue. Whether it is crys- 
talline alumina (corundum; AI2O3), a s first proposed by 
Glaccuml (11994ft. o r spinel (MgAl 2 4 ; IPosch et al.lll999l: 
Fabian et al.ll200lft , the groups favoring both candidates 
agree that an Al-0 stretch ing mode produces the feature 
(e.g.. lLebzelter et al.ll2006l ). 

Concentrating on this fact, we can infer that the pres- 
ence of the 13-/im feature indicates the presence of ele- 
mental aluminum. The 14 spectra with probable 13-/im 
features are distributed evenly over all metallicities down 
to [Fe/H]=— 0.71. Below this metallicity, all but two of 
the stars are naked, limiting any possible conclusions for 
that portion of the sample. The presence of aluminum 
at all metallici ties where it can be detected suports the 
conclusion in § 14.41 The absence of alumina in the Mag- 
ellanic samples was not a result of different abundances. 
We see alumina-based dust species in the globular sample 
at any metallicity where we see oxygen-rich dust. 

The difference between corundum vs. spinel as the 
carrier of the 13-/xm features boils down to whether 
or not simple oxide grains like MgO (or perhaps FeO) 
are actively bonding with the alumina. This differ- 
ence is impo r tant i n explaining the 20- /im features, since 
ISloan et ail ([2003). who favored corundum, suspected 
that the latter feature arises from crystalline silicates, 
while IPosch et al. (2002) have proposed an origin in sim- 
ple oxides ([Mg,Fe]0). In the former case, the presence 
of the 13- and 20-/iin features may point to the pres- 
ence of crystalline analogues to the amorphous alumina 
and silicates known to form in these dust shells, while 
in the latter, the two features may indicate a modified 
chemistry with enhanced oxides. Thus, a resolution of 
the question would help us understand the pathways fol- 
lowed by the grain chemistry in these dust shells. 

The presence of three spectra classified as "SXt" bol- 
sters the case for crystallinity as the origin of the 13- 
and 20-/tm features, since we can now see an accom- 
panying enhancement in crystalline olivine at 10 /im. 
These spectra are particularly intriguing, because they 
help to clarify the origin of the "three-component spec- 
tra", where narrow features at 11 and 13 /im are su- 
peri mposed on the broader 10 /tm silicate emission fea- 
ture ()Little-Marenin fc Little! I1988D . Here, we can see 
t he three componen ts more clearly than before. 

ISloan et al.l (120061). fo llowing theoretical arguments by 
Gail & Sedlmayr (1998), suggested that dust forming at 
lower densities may have a higher degree of crystallinity. 
Each photon absorbed by a dust grain will nudge recently 
accumulated atoms across the surface. If these atoms 
can be nudged enough times before the accumulation of 
an overlaying layer locks them into position, then they 
might fall into lowest energy levels represented by the 
lattice structure. 

A separate mechanism produces crystalline grains at 
higher mass-loss rates. The higher densities in the 
dust-formation zone push the condensation temperature 
higher, and the grains can anneal into a crysta lline struc- 
ture before they cool (see lFabian et al.l 120001 and refer- 
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Fig. 16. — The dust excess, after fitting and subtracting a stellar 
continuum, from the three spectra classified as "SY" due to their 
unusual 11-12 /an emission features. The vertical lines appear at 
wavelengths of 10, 11, 12, 13, and 20 (mi. 

ences therein). Thus, crystalline grains could form in 
cases of either particularly high or particularly low mass- 
loss rates. 

If indeed the appearance of features at 11, 13, 20, and 
28 /im is due to enhanced crystallinity, one could argue 
that this is evidence for a disk, which would retain the 
grains close to the star long enough for them to be ra- 
diatively annealed. We believe this is unlikely for most 
of the sources in our sample because disks would have 
higher optical depths and more cool dust than observed 
in these spectra. Generally, self- annealing during forma- 
tion at low densities is a better explanation than disks. 

Palomar 6 VI may be a different matter. If it is a clus- 
ter member, its bolomctric magnitude and period are 
inconsistent with the evolutionary track defined by the 
AGB stars in this sample (Figure [TOl and § I3.3[) . Its spec- 
trum is probably not produced by outflows from a normal 
AGB star. The presence of silicate absorption requires 
high optical depth, which could indicate the presence of 
an optically thick disk. Grains in such a disk would re- 
main in the star's vicinity longer than in an outflow, giv- 
ing them time to be annealed, which could explain the 
presence of emission features from crystalline silicates at 
23, 28, and 33 /tm. 

4.6. Unusual 11-12-^im Emission Features 

Figure [T|5] plots the three spectra classified as "SY" , 
due to a previously unrecognized dust emission feature 
peaking in the vicinity of 11-12 /tm. Two of the three 
sources, NGC 6553 V4 and Terzan 5 V9, also show 13- 
/tm features along with a feature at 20 /tm. Both spectra 
also show peaks in the emission at ~11.5 and ~12.2 /tm. 
Terzan 12 VI, while having general similarities to the 
other two spectra, differs in detail, with a single peak in 
the dust emission at ~11.1 /im and a relatively normal 
accompanying 18 /im feature. 
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Fig. 17. — Carbon dioxide emission bands in the spectra of three 
individual sources and the coaddition of six more. A spline has 
been fitted and removed from each spectrum, and the residual 
normalized to the maximum emission or absorption. The hor- 
izontal lines mark zero excess flux in each spectrum, and the 
vertical lines are at wavelengths of 13.48, 13.87, 14.98, 16.18, 
and 16.80 fim. The coadded spectrum at the bottom includes 
IC 1276 VI, NGC 5139 V42, NGC 5927 VI, NGC 6171 VI, 
NGC 6637 V4, and NGC 6760 V3. 

While the shape of the 11-1 2- /im emission feature is 
unusual, the flux ratios at 10, 11, and 12 /im place the 
three sources close to the silicate dust sequence in Fig- 
ure[T2](in the SE1-2 range, with Fu/F 12 < 1.05). Nor- 
mally, one would expect amorphous alumina to produce 
spectra with these flux ratios, but we have been unable 
to duplicate these spectra with optically thin dust shells 
consisting of only amorphous alumina, or even combi- 
nations of amorphous alumina and amorphous silicates. 
Optically thick models remain untested. It may be worth 
recalling that one of the three sources, Terzan 5 V9, is 
something of an enigma ( Figure [TO]), which could possi- 
bly be explained by the presence of a disk. 

4.7. Narrow Emission Features 

Many of the spectra show emission features at 14 /urn, 
but the majority of these fea tures are narrower than the 
14-/im features reported bv iSloan et al.l (|2006l . 2008). 
Only two sources show 14-/im features with similar po- 
sitions and widths to those in these previous papers: 
NGC 6356 VI and Palomar 6 VI. They are classified 
as 2.SE8f and 3.SBfx, respectively. The "x" for Palo- 
mar 6 VI indicates the presence of emission from crys- 
talline silicates at 23, 28, and 33 /im. The 14-/im dust 
emission feature tends to be associated with strong amor- 
phous silicate emission, as with NGC 6356 VI, or with 
crystalline features, as with Palomar 6 VI. 

The narrow 14-/im features seen in some spectra pre- 
vent a thorough search for the slightly broader 14-^m 
feature detected in NGC 6356 VI and Palomar 6 VI. 



The narrow 14-^m features are centered close to 13.9 /im, 
and some of the most pronounced examples are accom- 
panied by a second emission feature at 16.2 /im as well 
as absorption or emission at 15.0 /im. Figure [17] plots 
some examples. 

This combination of spectral features allows us to 
positively identify the bands as emission from CO2. 
Uusttanont et alj (|1998h first identified these bands in 
spectra from AGB stars obtained with the SWS. The 
detection of the same bands with the low-resolution mod- 
ules of the IRS is a bit of a surprise, because its spec- 
tral resolution is almost an order of magnitude lower 
than the SWS. 12 C02 produces narrow emission bands 
at 13.87, 14.98, and 16.18 /im, with the 14. 98 /im band 
shifti ng into absorption in some cases (see iCami et al.1 
2000). Some spectra also show the emission band seen in 
SWS data at 1 3.48 /im, as well as an additional band at 
16.8 /im, which ISloan et aD (|2003h argued was also from 
C0 2 . 

The repeated appearance of features at the right wave- 
lengths in our sample is convincing, even if many of the 
individual features in individual spectra remain in doubt. 
Given the noisy nature of the features, we have not at- 
tempted a quantitative analysis, though it is worth not- 
ing, that most of the sources included in FigurefTTlare not 
13-/im sources. If this result held up with better spec- 
tra of globular cluster variables in t he future, it would 
contradict the correlation found by Slo an et al.l (|2003l ) 
for oxygen-rich AGB variables in the Galaxy, where CO2 
emission appeared in spectra also showing 13-/xm fea- 
tures. 

4.8. Molecular absorption features 

Many of the spectra in the globular sample show clear 
bands from molecular gas shortward of 10 /im. The 
molecules include SO2 and H2O. 

Figure IT51 presents six IRS spectra showing absorption 
from SO2 at 7.3-7.5 /im, along with an SWS spectrum of 
UX Cyg, a Galactic AGB st ar and one of the first s ources 
in which SO2 was detected (jYamamura et al1 l2006) . SO2 
is the only sulphur-bearing molecule that has been de- 
tected in oxygen-rich AGB stars. One would not expect 
sulpher to be produced in low-mass AGB stars, and it 
would follow that the strength of the SO2 band should 
depend on the initial metallicity of the star. Thus, the 
presence of this band in metal-poor clusters like IC 1276 
and NGC 6441, with [Fe/H] = -0.6 or less, is unex- 
pected. 

The seven IRS and SWS spectra in Figure [TH] also 
show water vapor absorption at 6.4-6.8 /jm. Figure H"9l 
presents four more globular spectra with water bands, 
but no SO2, along with HD 32832, a Galactic M giant 
also observed with the IRS. This comparison spectrum 
also shows strong SiO absorption, which is generally ab- 
sent or weak in the globular sample. 

Figures [TH] and [TH] both contain synthetic absorp- 
tion spe ctra based on plane-p arallel radiative transfer 
models (iMatsuura et alJ |2002af) and line lists from HI- 
TRAN (jRothman et all 120091 ) These models include 
only the main isotopes, and the updated line lists pro- 
vide improved fits to the H2O structure at ~6.3 /tin 
(|Barber et al.| [2006) . The excitation temperature of the 
molecules is 1800 K. The absorption in the model in Fig- 
ure [18] is due entirely to water vapor and SO2 . In Fig- 
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Fig. 18. — Six IRS spectra showing SO2 absorption bands at 7.3— 
7.5 £im, along with a simple model and a comparison spectrum of 
a Galactic AGB star, UX Cyg, obtained with the SWS on ISO. 
The spectra are plotted in Rayleigh-Jeans units (A 2 -F„ \ 4 F\) so 
that the Rayleigh-Jeans tail of the blackbody function would be a 
horizontal line. All of the spectra also show absorption from water 
vapor at 6.4-6.8 fim. The simple model includes absorption from 
only H2O and S0 2 . 

ure !191 it is from water vapor and SiO. 

Figure [30] presents five more IRS spectra showing an 
apparent absorption band with a peak opacity ^6.35 /xm. 
The figure also includes a model showing water vapor in 
emission, which duplicates reasonably well the structure 
in the observed spectra. Thus, the "absorption" is ac- 
tually continuum between emission bands from H2O to 
ei ther side. 

iMatsuura et al.1 (|2002a[ ) found that LPVs tend to show 
water vapor in emission near the maximum of their pul- 
sation cycle, due to detached layers of water vapor in the 
extended atmosphere. Our five emission sources all have 
short periods, limiting our certainty of the phase during 
the IRS observations. The distinguishing characteristic 
of the emission sources is their periods. All five sources in 
our sample with periods between 40 and 140 days show 
water-vapor emission. The variables with longer periods 
d o not. 

iTsuiil (|2001ft detected water vapor in K giants, which 
are warmer than one would expect if the water vapor 
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Fig. 19. — Four IRS spectra of globular cluster variables show- 
ing water vapor absorption at 6.4-6.8 fim, but no SO2 absorption 
at 7.3-7.5 fim, a simple model, and HD 32832, a naked Galactic 
comparison source with a spectral class of M4 III observed by the 
IRS. The spectra are plotted in Rayleigh-Jeans units as in Fig. 1181 
The simple model only includes H2O and SiO. 
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Fig. 20. — Five IRS spectra showing a possible absorption feat ure 
at ~6.3 5 fim. The spectra are in Rayleigh-Jeans units, as in Fig. 1181 
and ll9l The model shows water vapor in emission, and it matches 
the structure seen in the data reasonably well. 
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were in hydrostatic equilibrium. Two of the emission 
sources in our sample are Cepheids, which are even 
warmer, and we can conclude that the water vapor 
around them is probably not in equilibrium. 

There has been some unpublished speculation about 
a possible dust emission feature at 6 /im and a possi- 
ble correlation with other features, such as the 14-/im 
dust emission feature. This "6-/xm emission feature" ex- 
tends from 6.0 to 6.5 /mi with a peak ~6.2-6.3 /im. It 
is particularly noticeable in the spectra of NGC 6441 V2 
and Terzan 5 V9 in Figure 1181 However, the complex 
spectral structure in the 6-/tm region created by molecu- 
lar absorption makes it more likely that this "feature" is 
really just the continuum between molecular absorption 
bands. 

5. OBJECTS OF NOTE 
5.1. Lyngd 7 VI 

The spectrum of Lynga 7 VI clearly identifies it as 
a carbon star, and our analysis strongly supports its 
membership in the cluster Q3.4p . While carbon stars in 
Ga lactic globular clusters are ra re, they are not unheard 
of. iWallerstein fc Knappl (|1998l ) listed three examples in 
NGC 5139 (oj Cen) and noted that all three were CH 
stars, which have probably beco me carbon-rich thr ough 
mass transfer in a binary system. ICote et all ()1997f ) dis- 
covered a candidate in the cluster M14; it too is a CH 
star. Lynga 7 VI may well be a similar object, which 
would explain how it could be a member of an old glob- 
ular cluster and still be a carbon star. 

For the sake of completeness, Table [4] includes a mass- 
loss rate for Lynga 7 VI based on its [6.4] — [9.3] color, 
the relation of col or to dust mass-loss rate defined by 
ISloan et "all (|2008f ). and an assumed gas-to-dust ratio 
of 200. Because of its potentially unusual evolutionary 
path, this mass-loss rate might not be useful for inter- 
comparisons among samples with different metallicities. 

5.2. NGC 5139 V42 

IMcDonald et al.l ([2009) obtained an N-band spectrum 
of NGC 5139 V42 about four weeks before our IRS ob- 
servation. Like us, they observed little dust in the spec- 
tr um, while earlier mid-i nfrared photometry (referenced 
bv lMcDonald et aT1l2009[ ) showed a clear excess at 10 /xm. 
They concluded that this change in the measurements is 
likely real. This spectral variation emphasizes the tempo- 
ral nature of these objects and the importance of obtain- 
ing statistically significant samples to average the varia- 
tions out. 

5.3. NGC 362 V2 and VI 6 

Our classification of a star as "naked" means that it 
does not show an identifiable excess in the 8-14 /im re- 
gion, based on its c ontinuum level in the 5-8 /im region. 
iBover et al.l (|2009bD recently detected a photometric ex- 
cess in two of our "naked" sources, NGC 362 V2 and 
V16, based on photometry of the continuum at shorter 
wavelengths. This excess is featureless (at the resolution 
of the photometry) , leading them to identify amorphous 
carbon as a likely suspect. Because the formation of 
CO would have left only oxygen to condense into dust 
in this oxygen-rich environment, the pr esence of carbon- 
rich dust would be most unexpected. IMcDonald et al.l 
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Fig. 21. — The dependence of -F11/-F12 (dust composition) and 
the 13-fim feature strength as a function of pulsation period and 
bolometric magnitude of the central star. The flux ratio -F11/-F12 
corresponds directly to the SE index, as indicated in Panel b. 

using the spectra presented here, show that iron 
grains provide the best explanation for the featureless 
excesses observed in these and other spectra. 

5.4. Cepheid variables 

Our globular sample includes four Cepheids. All four 
are type II Cepheids, with periods between ^20 and 
50 da ys, taken from the sample by [Matsunag a et all 
(2006). They could be classified as W Vir gims or 
RV Tauri stars, but the separation between these two 
groups is ambig uous especially among the globular clus- 
ter objects (See iMatsunaga et al.ll2009T ). 

RV Tau s tars are ofte n embedded within circum stellar 
dust shells ([Jural Il986lh INook fc CardeTIil <fl989h previ- 
ously reported an infrared excess from one source in our 
sample, NGC 6626 V17. However, the current IRS spec- 
trum docs not show a dust excess around this star, or 
the other three Cepheids. The difference may arise from 
temporal variations. 

It is interesting that the two Cepheids with longer pe- 
riods, NGC 6626 V17 and NGC 6779 V6, both show 
water-vapor emission in their spectra. The Galactic RV 
Tauri star R Sc uti (period 147 days) a lso shows water- 
vapor emission (jMatsuura et al.ll2002bl ). 

6. EVOLUTION ON THE AGB 

Figure |TD] compares the bolometric magnitudes and 
pulsation periods of the long-period variables in our glob- 
ular cluster sample to t heoretical evolutiona r y trac ks for 
low- mass AGB stars bv I Vassiliadis & Wood! (fl993f ). The 
tracks start at the theoretical beginning of the thermally 
pulsing AGB, and they are roughly consistent with the 
data from our sample. All of the data appear to fol- 
low a mutually consistent evolutionary track, although 
that track does not follow the theoretical track for either 
l.O-M or 1.5-M stars. The two sources which fall to 
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the right of and below the track are Terzan 5 V9 and 
Palomar 6 VI, and their location away from the rest of 
the data has already led us to conclude that either they 
are not actually cluster members, or they are not normal 
mass-losing AGB stars. 

In their study of 47 Tuc, iLebzelter etaD (j2006| ) found 
that their IRS spectra were consistent with the onset of 
dust formation at a luminosity of ^2000 L , which corre- 
sponds to Mm ~ —3.5, followed by a switch in pulsation 
mode from an overtone to the fundamental. The dust 
spectra in their sample showed strong 13-/Ltm features 
at this stage, which then grew comparatively weaker as 
amorphous silicates begin to dominate. They empha- 
sized the association of a stronger 13- /im feature with the 
fainter and less evolved stars in their sample. Our study 
cannot address the question of a switch in pulsation mode 
because all of our long-period variables appear to be pul- 
sating in the fundamental mode, which is consistent with 
all of our data being brighter than Mboi ~ —3.5. 

Figure [2T1 tracks the dust properties in our sample with 
pulsation period and bolometric magnitude. The top 
panels show that the 13-/im feature is most pronounced 
for pulsators with periods around ^200 days and bolo- 
metric magnitudes from ~— 3.5 to —4.0. Thus, our glob- 
ul ar sample genera l ly con forms to the scenario described 
by ILebzelter et all ()2006l ) . with the 13-/Ltm feature grow- 
ing weaker and the contribution from amorphous silicates 
increasing as the star evolves along the AGB. We do not 
see any evidence for sudden shifts in spectral properties, 
but then our sample shows no shift from overtone to fun- 
damental mode, as seen in 47 Tuc. 

The bottom panels of Figure |2"T1 show little dependence 
of the corrected flux ratio F\\/Fn with either period or 
bolometric luminosity. We have noted previously that 
Galactic supergiants and the luminous AGB stars and 
supergiants observed in the Magellanic Clouds show little 
sign of amorphous alumina in their spectra. The globular 
sample, though, is restricted to lower luminosities, and 
within that range, amorphous alumina and amorphous 
silicates can dominate the shells with little apparent de- 
pendence on either period or luminosity. 

7. SUMMARY 

We have presented the infrared spectra of a sample of 
39 variable stars in 23 globular clusters. The sample in- 
cludes four Cepheid variables, none of which show dust in 
their spectra, and 35 long-period variables, of which we 
can confirm 31 as members of their clusters. These stars 
are either naked, show emission from oxygen-rich dust, or 
in one case, show emission from carbon-rich dust. The 
variety of dust species is impressive, and we see amor- 
phous silicates, amorphous alumina, emission from crys- 
talline grains at several wavelengths, and three spectra 
which we have not yet been able to characterize. The 
spectra are also rich in molecular features, including ab- 
sorption from H2O and SO2, and emission from H2O and 
C0 2 . 
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The main objective of this paper is to probe how the 
quantity and composition of the dust depend on metal- 
licity. Our sample shows dust emission over the range 
—0.97 < [Fe/H] < —0.08. Across this range, stars with 
lower metallicities generally have less dust in their cir- 
cumstellar shells. This trend is most readily appar- 
ent when plotting (dust) mass-loss rate as a function 
of the dust composition (as quantified by the flux ratio 
F11/F12), but it is also noticeable when plotting mass- 
loss rate vs. bolometric magnitude. The lack of a metal- 
licity trend vs. pulsation period may arise from depen- 
dencies of the period on metallicity. 

We find contributions from alumina-rich dust at every 
metallicity at which we see silicates. This result indi- 
cates that the lack of alumina in the Magellanic samples 
published thus far is most likely a selection effect, since 
those samples do not include objects as intrinsically faint 
as those published here, and the brighter objects tend to 
show mostly silicates. 

Three spectra have a 10 /im-emission feature split into 
components at 9.7 and 11.3 /im, along with a strong 13- 
(im feature. The splitting at 10 /zm is reminicent of two 
spectra in the LMC which [Sloan et all ([20081 ) fitted with 
crystalline analogues of the amorphous silicates which 
produce the broad 10- /urn feature. These spectra support 
the case that the 13- /im feature arises from the crystalline 
analogue of the amorphous alumina seen in several other 
spectra. 

This initial paper has only scratched the surface of a 
rich set of spectroscopic data. A careful mineralogical 
analysis of the dust features should help address some 
of the questions left unanswered here. In addition, the 
spectra are rich in molecular features, both in the 5-8-/xm 
region and in the 13.5-17 /jm region. Detailed modelling 
will give us insight about the behavior fo these molecules 
in metal-poor environments. 

While Spitzer has exhausted its cryogens and the IRS 
no longer operates, there is still a strong need for further 
infrared spectra of evolved stars in globular clusters. The 
sample here is only large enough to point to trends of 
fundamental importance. The need for more data and 
larger samples is clear. 
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